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ANGLE-OF-ARRIVAL  FLUCTUATION  IN  BEAM 
PROPAGATION  THROUGH  TURBULENT  ATMOSPHERE 

Song  Zhengfang 

Anhui  Research  Institute 
of  Fine  Optical  Machinery 
Chinese  Academy  of  Sciences 
Hefei  City,  P.O.  Box  1125, 

230031 

ABSTRACT.  In  this  paper,  starting  out  near-axis 
approximation  and  considering  the  correlation  of 
intensity  fluctuation,  one  derives  a  universal 
formula  (Eg,  14)  of  angle-of-arrival  (AOA)  fluctuation 
variance  under  the  assumptions  of  a  Markov  random 
process.  For  plane  and  spherical  wave  propagation, 
analytic  expressions  of  are  obtained.  Some  experimental 
data  are  used  to  compare  with  analytic  results .  It  shows 
that  both  are  identical . 

KEY  WORDS:  Turbulent  flow;  propagation  through  atmosphere, 
angle  of  arrival . 


1 .  Introduction 

After  propagation  through  the  atmosphere,  the  location  of 
optical  images  on  the  focal  plane  of  receiving  telescopes  shows 
random  trembling  due  to  perturbation  of  turbulent  flow.  This  is 
the  result  of  so-called  fluctuations  of  beeim  angle  of  arrival. 

In  the  past  two  decades,  a  huge  amount  of  theoretical  and 
experimental  research  was  focused  on  this  subject,  yet  most 
studies  were  devoted  to  the  characteristics  of  the  weak 
fluctuation  zone.  As  for  the  angle  of  arrival  in  the  strong 
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fluctuation  zone,  little  was  known  thus  far.  Experiment  shows 
that  the  weak  perturbation  theory[13  can  be  used  to  explain 
satisfactorily  the  experimental  laws  for  weak  turbulent  flows. 
However,  it  will  remarkably  deviate  from  the  experimental  values 
with  strong  turbulent  flows  or  long  distances[2] . 

It  is  well  known  that  under  strong  turbulent  flows,  the 
Markov  approximation  can  be  used  to  treat  wave  propagation 
equations  with  satisfactory  results.  However,  while  considering 
the  correlation  functions  of  intensity,  it  is  possible  that 
correct  results  still  cannot  be  obtained  if  the  second-order 
moment  of  fields  is  not  taken  into  calculations.  In 
reference[ 3] ,  we  tentatively  discussed  the  angle  of  arrival  by 
using  Markov  approximation  and  wave  complex  amplitude  normal 
distribution  approximation.  In  this  paper,  we  will  further 
discuss  the  expressions  of  angle-of -arrival  variance,  which  are 
applicable  to  the  entire  fluctuation  zone.  In  addition,  we  will 
also  detail  the  plane  wave  and  spherical  wave  propagation  and 
compare  the  existing  theories  with  experimental  results.  The 
comparison  suggests  that  our  practice  here  is  rational . 

2 .  Foirmula  Derivation 

The  basic  expression  for  the  angle  of  arrival  is[3] 

(1) 

where  Pq  is  total  flux  received,  L  is  propagation  distance, 
vp  is  the  gradient  operator  of  horizontal  vector  p  ,  ni  is  the 
fluctuation  value  of  refractivity,  =  u(^,p)  is 

the  complex  amplitude  of  waves,  and7’(p)is  the  transfer  function 
of  the  receiving  telescope.  From  Eq.  (1),  the  expression  for  the 
angle-of-arrival  variance  can  be  derived  as  follows: 
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(2) 

(3) 


where  <...>  is  the  system  integration  average  selected.  Based  on 
Markov  approximation  and  the  assumption  that  the  fluctuation  of 
the  angle  of  arrival  mainly  comes  from  phase  fluctuation,  the 
following  can  be  verified; 
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2k^Pl 


Pid-PiVl^si^  ,  p)<7(  ^,p,m  .ft)>  T(fi^)T(pd 


(4) 


where  k=2it/X;  X  is  wavelength,  M^,'p)  is  the  phase  structure 
function,  which  already  contains  the  integration  of  distance. 


To  conduct  a  general  investigation,  a  wave  field  in  Gaussian 


form  is  adopted.  At  this  instant,  u(^,^)  =  i/„exp(- 


2cr 


t) 


o-=-a,-’-i± 

F 


(5) 


where  Oo  is  the  equivalent  radius  of  laser  beam  at  the 
transmitting  end,  and  F  is  the  curvature  radius  of  the  wavefront. 
Next,  the  transfer  function  T(  p  )  is  also  believed  to  be  in 
Gaussian  form 


T(p)  =  Toexp(-  -^> 
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where  To  is  penetration  coefficient,  otx  is  the  radius  of  the 
receiving  telescope. 


Under  the  above  condition,  it  is  also  necessairy  to  know  the 
form  of  the  intensity  correlation  function  </(^,Pi)^(s.P:)>  ii' 

solving  Eq.  (4),  which  used  to  be  selected  as 

This  is  approximately  effective  in  the  weak  fluctuation  zone,  but 
in  the  entire  fluctuation  zone,  the  correlation  of  intensity 
fluctuation  must  be  taken  into  account.  We  confirmed[4]  that  it 
would  be  appropriate  to  adopt  the  following  formula: 

</(4,pi)/(^P2)>  =</(^.pO></(^,P2)>+ir,(^p,,p,)i^  (7) 

where  PM,p\,'p^  second-order  moment  (or  coherent 
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function)  of  the  field,  whose  expression  was  given  by  Wang  and 
Plonus[5]  by  using  the  quadratic  approximation.  Here,  the 
coherent  light  source  can  be  rewritten  as  follows: 


.  Pc .  )-exp  ( -  ibp,  •  ft - ^  -  dfi) 

«  can 


K-  Y  (^\  +  Pd 


P*=P2-P} 


+  4-  +(-4-)(i-4)' 


C=(I-  -^  )'+(]+£,  V' 

*  # 

)  +  (4- )’  +  (!+ 


loo 

fif= 


where  f  is  the  Feiner  number  of  the  transmission  field,  g  is  the 
specific  value  between  beam  width  and  the  coherent  length  of 
spherical  wave  po  or  the  relative  beam  width  in  the  turbulent 
atmosphere.  The  expression  of  po  is 

Po=(0.546Clk-L)-^>^  (9) 

where  Cn  is  the  refractivity-f luctuation  structure  constant, 
which  indicates  the  intensity  of  turbulent  flows. 


In  f  ^=0  taken  as  the  average  intensity 

<  /(s,ft)  By  replacing  Ri=  Pi  and  R2=  Pi  in  Eq.  (8)  and 

introducing  variables  R=Y(ft+py  and  P=ft-ft  /  then 

inserting  all  these  together  with  Eq.  (6)  in  Eq.  (4),  the 

following  can  be  derived: 

2_  “o  /  /  ^<T’2  ill 


(10) 


.  -MV)) 

In  the  above  equation,  we  already  let 

ar-=(ca^"'  +  af2 

d'  =  d  +  (4ci^j.)-' 

By  integrating  variable  R  in  Eq.  (10)  first  and  then 
calculating  Po,  Eq.  (10)  will  change  to 

f*  0^ 

a;=(2k-)"'(a7  +  a^"'  dppV;jg)s(<‘ ,p)  •{  exp(-  )  +  exp(“  Yr 

Jo  *  “ 


(11) 
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In  the  calculations,  we  already  allowed  the  following: 

Like  ai,  02  also  has  a  dimension  of  length.  It  can  be  believed 
that  they  have  equivalent  integrating  radii  of  different  values. 
Their  difference  is  that  aj  not  only  is  related  to  beeun  radius 
and  receiving  radius,  but  also  is  determined  by  the  intensity  of 
turbulent  flows. 


Eq.  (11)  is  a  general  expression  of  Gaussian  beam  angle-of- 
arrival  variance.  Since  this  kind  of  beam  has  a  very  complicated 
phase  structure  function,  there  is  generally  no  analytic  solution 
at  all.  Therefore,  to  understand  the  characteristics  of  of , 
appropriate  measures  have  to  be  taken  in  response  to  specified 
conditions.  As  an  example,  the  following  discussion  is  devoted 
to  plane  wave  and  spherical  wave  propagation  so  as  to  ascertain 
the  appropriateness  of  our  measures  here. 

3 .  Plane  Wave  Propagation 

Both  theoretical  analysis  and  practical  examination  suggest 
that  the  phase  structure  function  derived  based  on  the  classical 
perturbation  method[l]  is  also  applicable  to  the  strong 
fluctuation  zone,  which  serves  as  the  foundation  for  our 
discussion . 


The  phase  structure  function  of  plane  waves  can  be  expressed 


in  the  following  equation: 

mL  p  ® 

JO 


,p)  =  4K^k^{dA 

Jc  Jo 


k 

oT  Sin 


■) 


(12) 

where  $n(K)  is  the  density  function  of  turbulence  spectrum,  Jo(x) 
is  zero-order  Bessel  function.  Considering  ~  Jo(/>A)]  =  A-Jo(Ap)  ^ 

in  Kolomogorov  turbulence  modeling.  Thus,  Eq.  (11)  changes  to 
<I>,(K)  =  0.033C;A:""^’ 


^2_  0.0667r 


diCUo 


dpp[exp{-  ^)  +  exp(-  )] 
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dKK-^Jo{Kp){\  +  sin  -^  ) 


JO  (13) 

The  first  part  of  integration  of  variable  K  is  easy  to  make, 
while  its  second  part  can  also  be  solved  through  the  ready 
integration  formula  after  x=K^$/k  is  replaced.  Based  on  this 
outcome,  the  integrating  operation  can  be  done  for  p,  resulting 
in 

ff2=2.564(2a, 


’P 


1  +m 


W 


•Iin(i-’^[>,(-  -j  ,  1 ;  1;  «0)  +  m  2f,(-  ,  1 ;  1  ;m»f2)l)} 


(14) 


where  we  already  let 


m= 


Q= 


ka] 


Im( • )  suggests  selecting  the  imaginary  part  of  the  function  in 
the  parentheses,  2Fi(ai/a2/t>/X)  is  the  hypergeometric  function, 
Q  is  the  Feiner  niomber  of  the  integration  field. 


Eq.  (14)  is  the  basic  equation  for  the  angle-of-arrival 
variance  of  plane  waves,  which  can  be  solved  as  long  as  the 
relationship  between  turbulence  intensity  and  distance  is 
ascertained.  To  save  space,  here  we  only  discuss  the  uniform 

In  this  case,Cj(*)is  a  constant,  so  Eq.  (14)  can 


turbulent  flows, 
be  rewritten  as: 


«7;=2.564C„^L(2a,)-'^(-^j^)|l  + 


(15) 


As  far  as  the  collimated  laser  becim  with  limited  space  is 
concerned,  it  can  satisfy  the  condition  PoXXo  when  turbulent 
flows  are  extremely  weak.  In  addition,  in  the  far  field, 

(/c/s«  1)  ’,  and  the  specific  value  of  the  integrating 

radius  m=l .  Meanwhile,  the  condition  of  the  Feiner  number  of  the 
integration  field  Q<1  can  also  be  established.  Therefore, 

Eq.  (15)  can  be  simplified  as: 

(16) 

ol=2.56ClL(2(Xry''^ 
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agrsGs  closoly  with  th©  jrssults  of  th©  w©&k  p©]rt'uirb3,tion 
th©ory[ 1] . 


For  th©  conv©n.i©nc©  of  calculations,  th©  hyp©rg©om©tric 
function  can  b©  ©xpr©ss©d  in  s©ri©s.  Th©n  Eq.  (15)  will  b©com©  a 
simpl©  formula  as  follows:  di=(^{Ci) 

.  1.160''^  (17) 


/(«)= 


1  +m 


o  +  i 


+  m 

00 


,5/6 


•[  +m^*')+ 


»»0 


l!“0 


(2fi-17/6)(2n-ll/6) 

- 2n(2n-l  j  ’  ’ 


(2n-ll/6)(2n-5/6) 

2n(2n+l) 

6o= -0.2233 


Fig.  1.  f(J2)  as  a  function  of 

Fig.  1  lists  part  of  calculations  conducted  for  function  f(fl). 

It  is  known  from  th©s©  calculations  that  in  general  situations, 
al  is  related  not  only  to  turbulence  intensity  Cr,  ,  propagation 
distance  L,  and  receiving  aperture  otj  as  usual,  but  is  also 
related  to  the  Feiner  parameters  of  the  receiving  field  and  the 
specific  value  of  the  integrating  radius.  This  is  probably 
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because  the  effect  of  turbulent  flows  on  light  beams  depends  on 
the  relative  relations  between  beam  scale  and  turbulence  scale. 
While  the  weak  perturbation  theory  deals  with  infinite  plane 
waves,  whose  scale  is  much  larger  than  turbulence  scale,  and 
therefore,  it  is  not  involved  in  such  relations  as  mentioned 
here . 


4 .  Spherical  Wave  Propagation 


The  phase  structure  function  of  spherical  waves  is 


-^o(  )lcosH 

Jo  _ . 


] 


(18) 


Noting  that  j ,  _  ^  2^  j,(-^  )  and 

cos^x=  —  (1  +cos2;c)  /  under  the  condition  of  the 

Kolmogorov  turbulence  spectrum,  Eq.  (11)  changes  to 


a:= 


0.066;r 

L- 


j^^PP{exp(-  -^)  +  exp(- 


OT  +  t 
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{ (JP)-IP+  lOM  ( — ^  .  .^'/op^r.-i/o  C  ,  _L  ■  1  ■ 
2^r(5/6)  ^  L  ^  2 


(19) 


By  further  integrating  variable  p  and  considering  that  Qfci.>>ao  can 
always  be  satisfied  in  the  case  of  spherical  waves,  i.e.,  ai=Qtx, 
the  general  expression  of  the  angle-of-arrival  variation  of 
spherical  waves  is  as  follows: 


fff=2.564L~*^ 


iCiji 

—  )  +  I-l- - ^ 


1/6 


•  Re(f-'‘[  :F,(  ,  1 : 1 ;  jF,(  ^  ,  I ;  1 ;  )])} 

where  we  already  let  ^ 


(20) 


/Wt*  ~ 


^  2L 


Just  in  the  case  of  plane  waves,  Eq.  (20)  can  also  be  solved 
as  long  as  the  function  form  of  C^((J)is  known.  For  the  uniform 
optical  path,  we  obtained 

= 0.96 1 5CjL(2a7.)~ '  V(nT) 


/(fir)= j  1+0.09477  ^ 


•[Gjj(«i2r|-  -j  ■  0;  -|.0,0)  +  m^Gi:5(2m^Q^|-  ,  A  ,0;  ,0,0)]j| 


(22) 


where  G”^{x\a,;b,)  is  the  Meijer  G  function.  The  relationships 
between  the  G  function  and  the  hypergeometric  function[7]  were 
already  made  use  of  in  solving  Eq.  (21). 


The  G  function  is  extremely  complex  in  nature.  However, 
with  it  we  can  analyze  the  properties  of  function  f(QT)  in  the 
near  field  (Qt>1)  and  far  field  (Qt<i)*  Under  the  conditions  of 
the  near  field,  mx^lt  /  while 


(^>1) 


Then,  Eq.  (22)  can  be  rewritten  as 

Aar)=  |l  +0.79369(1  +m?^-'Re^-''’[  J,(l .  ±  -  2  ;  jQ;') 

+  mr,/.(l.  -  y  ;<(m,flrr')l-0.74698iQV(  ^  ;-U.  ;jO;>) 

Under  the  conditions  of  the  far  field,  mT®*  p^ACI^t  ,  while 


Gij(x|- 


n 


;x) 


+  r(4  )r( 


-^)r(-g-):F,(  — 


^■±X) 
’  6  ’  6 


(x<l) 


Thus,  Eq.  (22)  becomes 


/(nr)= 


\+mf 

*  +mj 


{ 


1+  1.0269 


( 


+  '«r3F,(Jl  ^ 

+  "«r2^i(  1  j  ^  (Qr<  1) 


(23') 


Once  the  Feiner  number  of  the  receiving  field  and  the 
coherent  length  of  spherical  waves  Po  are  determined,  function 
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can  be  derived.  Fig.  2  and  Fig.  3  present  part  of  calculations 
of /(^t)  .  In  Fig.  2,  under  the  limitation  of 

£V»1  ,  because  at  this  point,  ntr  is  equal  to  QV  .  In  Fig. 

3,  only  a  few  numbers  of  mT<l  are  selected,  because  in  the  far 
field,  the  condition  of  Po<aj~  can  always  be  satisfied. 


Fig.  2.  f(QT)  as  a  function  of  Qy  Fig.  3  fccr)  as  a  function  of  Siy 
(QT>1)  (QT<1) 

To  demonstrate  the  effectiveness  of  the  above  treatment,  we 
made  a  comparison  between  the  experimental  results  listed  in 
reference[2]  and  the  calculations.  The  conditions  of  the 
experiment  are:  L=1750m,  A,=0.5|im,  06^=2. 36cm,  point  light  source. 
According  to  the  definitions  given  in  this  paper,  we  obtained 
flT=2;  we  also  obtained  f(QT)  =0-804  from  Eq.  (23).  Fig.  4  shows  a 
comparison  between  the  experiment  and  the  theory,  where  the 
straight  line  1  is  l.03C„L''^(2ctr)‘''‘’  /  derived  from  the  weak 

perturbation  theory,  while  the  straight  line  2  is  the  result  of 
this  paper.  This  figure  indicates  that  the  expected  value  in  the 
weak  perturbation  theory  already  deviates  greatly  from  the 
experimental  value  at  a>20^rad,  while  the  results  given  in  this 
paper  conform  to  the  actual  situation  to  a  greater  degree.  It  is 
noteworthy  that  when  Oa;,20ti.rad /  equivalent  to 

C-  =( 2~6 )xlO'^%'^^^,  being  a  mid  to  strong  turbulence  state.  It 
can  be  further  expected  that  Eq.  (23)  is  possibly  applicable  to 
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the  strong  fluctuation  zone. 


5.  Conclusions 

By  using  the  Markov  approximation,  and  the  quadratic 
approximation  of  average  intensity,  we  acquired  the  general 
expression  of  the  angle-of -arrival  variance  and  derived  the 
analytic  formula  of  for  plane  waves  and  spherical  waves.  The 
practicability  of  these  results  is  surely  limited  due  to  the 
foregoing  approximation. 


Fig.  4.  Comparison  of  theory  with  experiment 
KEY:  1  -  small-perturbation  theory  2  -  our 
calculation 

However,  the  Markov  approximation  is  a  traditional  method  dealing 
with  wave  propagation  through  turbulent  media  and  has  been 
supported  by  experiments [ 1 ] .  As  for  the  square  approximation  of 
average  intensity,  it  serves  as  an  ideal  approximation  according 
to  Leader[8].  As  a  matter  of  fact,  our  calculation  results  of 
plane  waves  proved  to  conform  to  the  mature  theory  with  regard  to 
weak  turbulent  flows,  and  our  calculation  results  of  spherical 
waves  were  supported  by  experiments,  which  suggests  that  the 
above-mentioned  assumption  is  reasonable. 
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In  summary,  by  comparing  our  calculations  with  the  mature 
theory  and  the  experiment,  it  is  demonstrated  that  the  results 
obtained  in  this  paper  are  rational  and  possibly  are  applicable 
to  the  entire  fluctuation  zone.  Despite  this,  our  results  are 
expected  to  be  further  verified  due  to  insufficient  experimental 
data  and  particularly  those  under  different  conditions. 
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